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Abstract: A key insight used in developing the theory of Causal Dynamical Trian- 
gulations (CDTs) is to use the causal (or light-cone) structure of Lorentzian manifolds 
to restrict the class of geometries appearing in the Quantum Gravity (QG) path in- 
tegral. By exploiting this structure the models developed in CDTs differ from the 
analogous models developed in the Euclidean domain, models of (Euclidean) Dynam- 
ical Triangulations (DT), and the corresponding Lorentzian results are in many ways 
more "physical". 

In this paper we use this insight to formulate a Lorentzian signature model that 
is analogous to the Quantum Regge Calculus (QRC) approach to Euclidean Quantum 
Gravity. We exploit another crucial fact about the structure of Lorentzian manifolds, 
namely that certain simplices are not constrained by the triangle inequalities present in 
Euclidean signature. We show that this model is not related to QRC by a naive Wick 
rotation; this serves as another demonstration that the sum over Lorentzian geometries 
is not simply related to the sum over Euclidean geometries. By removing the triangle 
inequality constraints, there is more freedom to perform analytical calculations, and in 
addition numerical simulations are more computationally efficient. 

We first formulate the model in 1 + 1 dimensions, and derive scaling relations for 
the pure gravity path integral on the torus using two different measures. It appears 
relatively easy to generate "large" universes, both in spatial and temporal extent. In 
addition, loop-to-loop amplitudes are discussed, and a transfer matrix is derived. We 
then also discuss the model in higher dimensions. 
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1 Introduction 



Prior to the formulation of Causal Dynamical Triangulations (CDTs) in reference [1], 
and further developed in [2-6], the two main approaches to simplicial Quantum Gravity 
(QG) were those of (Euclidean) Dynamical Triangulations (DT) and the (Euclidean) 
Quantum Regge Calculus (QRC), also known as (Euclidean) Fixed Triangulations 
(FT). Both of these approaches are formulated in the Euclidean sector of QG, that 
is, they are sums over geometries which have Euclidean signature. In two dimensions, 
it is possible to completely solve the DT path integral by using matrix model tech- 
niques (see [7] and references within, and [8] for the generalization to CDT), however 
the results are completely unphysical: There exist only two phases of the quantized 
geometry, a crumpled phase and a polymerized phase [9] ; this rather strongly suggests 
that DT does not have a physically relevant classical limit. This problem was remedied 
in CDTs by shifting to a sum over Lorentzian-signature geometries which exhibit a 
"causal" structure. That is, one sums only over that subset of Euclidean geometries 
that is compatible with the existence of a Lorentzian signature metric. The key point is 
that CDTs exploit the causal structure present in Lorentzian geometry, and in so doing 
lead to a more physically realistic model of Quantum Gravity. (For related comments 
in continuum quantum gravity, and further background references, see [10].) 

In counterpoint, the QRC approach to Euclidean QG is rather hard to work with 
analytically [11] (except in the weak field approximation), owing to the fact that each 
simplex in the triangulation must satisfy generalized triangle inequalities. However, in 
the Lorentzian domain this difficulty disappears — completely so in 1+1 dimensions and 
largely so in higher dimensions — provided that the simplices in the triangulation are 
chosen such that there are both time-like and space-like edges. QRC in the Lorentzian 
domain has been discussed before in reference [12], where the inverse free propagator 
was calculated, however the disappearance of certain triangle inequality constraints was 
not noted or discussed. 

In section 2, we demonstrate that for Lorentzian triangles with one space-like and 
two time-like edges (or vice- versa), the relative magnitudes of the edge lengths are 
completely unconstrained. The Lorentzian triangle inequalities are, in this particular 
context, vacuous. Exploiting this fact, in section 3 we formulate a model of Lorentzian 
quantum gravity in 1+1 dimensions and demonstrate that it is not simply related to 
the QRC model in 1+1 dimensions by Wick rotation. In section 4 we derive scaling 
relations for the pure gravity model on the torus using two different measures, one of 
them corresponding to the DeWitt measure in 2 dimensions. In section 5 we discuss 
loop-to-loop amplitudes and derive a transfer matrix for the model. In section 6 we 
will discuss the non-occurrence of spikes in this model. In section 7 we discuss the 
configuration space of the theory in higher dimensions. We then conclude with a brief 
discussion. 
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2 Lorentzian Triangles 



In formulating the Lorentzian-signature Fixed- Topology Triangulation (LFT) model in 
1+1 dimensions we consider Lorentzian triangles which have one space-like edge and 
two time-like edges. Unlike triangles in Euclidean space, these triangles (as well as ones 
with two space-like and one time-like edge) do not have their edge lengths constrained 
by inequalities. Without loss of generality consider a Lorentzian triangle with one 
space-like edge length S, and two time-like edges of lengths T\ and T 2 in the coordinate 
system shown in figure 1. The third vertex of the triangle lies at the point (t, x), which 
must be in the future light-cone of both of the other vertices, i.e. t > and \t\ > \x\ 
(the case where it lies in the past light cone follows in an identical manner). Therefore 
we have 

T 2 = t 2 -x 2 + S{t-x) >0, (2.1) 
T 2 = t 2 -x 2 + S(t + x) > 0. (2.2) 



(*,*) 
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Figure 1. A Lorentzian triangle with space- like edge length S, and time-like edge lengths 
Ti, T 2 . 



These relations can be inverted to yield 

t = ^ (VW + [Ti - T 2 ] 2 ){S* + [T\ + W) 



X = ^( T 2-Ti)- 



(2.3) 
(2.4) 
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Both equations (2.3) and (2.4) have a real solution for any value of S G (0, oo) and Ti, 
T 2 G [0, oo); thus the edge lengths are completely unconstrained. For given values of S, 
Ti, and T 2 , we note that in the coordinate system under consideration the third vertex 
is located at the intersection of the hyperbolae given in equations (2.1) and (2.2). Using 
equation (2.3) we can easily find the area of a Lorentzian triangle in terms of the edge 
lengths S, Ti and T 2 : 

A L (S } T 1} T 2 ) = 1 -xSx f| +tj, 

= -^{S 2 + [Ti - T 2 ] 2 )(S 2 + [Ti + T 2 ] 2 ) , (2.5) 
= ^S 4 + T? + Ti + 2S 2 T 2 + 2S 2 T 2 - 2T 2 T 2 . (2.6) 

Equation (2.6) is closely related to, but not identical to, Heron's formula which describes 
the area of a Euclidean triangle in terms of its edge lengths a, b, c by: 

A E (a, b, c) = ^V-a 4 - 6 4 - c 4 + 2a 2 6 2 + 2a 2 c 2 + 26 2 c 2 . (2.7) 

The relation between these quantities is what one would expect by a simple Wick 
rotation, t — > it: 

A L {S,T X ,T 2 ) = —iA E (S, iT-i, iT 2 ), (2.8) 

where the factor of —i is due to the fact that we should be multiplying Lorentzian area 
by a factor of a/ det(r/) = i. Again, it is clear that these Lorentzian triangles are not 
constrained by inequalities, for while in equation (2.7) the area becomes pure imaginary 
if the Euclidean triangle inequalities are violated, in equation (2.6) the area is real for 
any choice of S, Ti, T 2 . 

Another quantity that is of interest in studying simplicial gravity is the angle 
between two edge lengths. For the angle between two space-like edge lengths these are 
found using the standard trigonometric relations, however finding the angle between 
a space-like edge and a time-like edge or two time-like edges requires more care. For 
configurations such as in figure 1 these angles will be complex numbers (boosts). A 
consistent way of assigning angles to these triangles is described in reference [13]. We 
denote the angle between S and Ti to be 61, between S and T 2 to be 6* 2 , and the angle 
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between T% and T 2 to be tp. These angles are given by: 

^A L {S,T X ,T 2 ), (2.9) 

^rA L (S,T lt T a ), (2.10) 

£=- A L (S,T U T 2 ). (2.11) 

L XJ-2 

There are few things to notice here: First is that cos 9\ and cos 9 2 are purely imaginary, 
and therefore = 9? [#2] — 7r /2, and second is that cosy? > 1 so that §ft[</?] = 0. 

Equipped with these notions of Lorentzian-signature triangles, we shall now define 
a theory of Lorentzian-signature Quantum Gravity in 1+1 dimensions, which will be 
closely analogous to QRC However as we will soon demonstrate, these theories are not 
related by any simple Wick Rotation. 
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3 1+1 Lorentzian Fixed Triangulations 

Classically two-dimensional Einstein gravity is trivial, since the action is a topological 
invariant, the Euler characteristic, and is therefore constant with respect to local de- 
grees of freedom in the theory. This is equivalent to the fact that in two dimensions 
the Einstein tensor vanishes identically. However two-dimensional QG is non-trivial, 
and due to the topological nature of the Einstein-Hilbert term can often be solved an- 
alytically. Two dimensional simplicial QG was the framework within which CDT was 
first formulated [1], and has also been investigated using QRC [14-16]. Here we will 
consider pure two dimensional QG, that is, the path integral: 

Z = J V[g] e iS ^ , S[g) = A J d 2 x^g . (3.1) 

In the Regge Calculus [17, 18], a discrete approach is taken in evaluating equations 
(3.1). The manifold is taken to be piecewise flat and is composed, in two dimensions, 
of triangles glued along shared edges. The metric degrees of freedom are described 
by the edge-lengths of each triangle, and curvature is localized at the vertices of the 
triangles. We shall consider a model in which the manifold under consideration either 
has cylindrical topology S 1 x R, where the S 1 are space-like surfaces parameterized by 
"time", or the toroidal topology S 1 x S 1 . 

3.1 Triangulations 

This manifold has a simple triangulation the "traditional" version of which (in analogy 
with common conventions in Euclidean signature) is shown in figure 2. Here the hori- 
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Figure 2. An 8x8 triangulation of the cylinder and a vertex in the triangulation with its 
associated edge lengths. (Traditional version.) 

zontal edges are taken to be space-like, and the diagonal and vertical lines are taken to 
be time-like. For the cylinder the time-like boundaries are identified, and for the torus 
the space-like boundaries are also identified. A physically equivalent but more "sym- 
metric" version of the triangulation is presented in figure 3. The simplices in either of 




Figure 3. Regular triangular lattice. (Symmetric version in terms of S, T\ and T2.) For 
each triangle one edge is spacelike, two are timelike (or at worst null). 

these triangulations are the Lorentzian-signature triangles described in section 2. We 
take there to be N s edge-lengths in the spatial direction, and N t in the time direction, 
for a total of A^ = 2N s N t triangles in the spacetime. Note that by adopting a fixed- 
topology triangulation of this type one is automatically excluding, at the kinematical 
level, the possibility of a polymerized phase. By this we mean that the polymer-like 
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phase found in DT [9] is not realizable in this model, however this does not exclude 
the possibility of a "rough" phase, (that is, a "fractal" or "crumpled" phase), similar 
to that found in QRC [14-16]. Before we rewrite equations (3.1) in their simplicial 
form we point out a few differences between two dimensional Euclidean and Lorentzian 
Regge calculus. 

3.2 Metric: traditional 

For Euclidean triangulations, the flat metric inside a triangle with edge lengths 1%, I2, 
and Z3 (in a coordinate system where l\ is the unit vector in the 1 direction and I2 is 
the unit vector in the 2 direction) is given by: 



Note that this satisfies det(g) = 4A| ; (Z 1 , Z 2 , Z 3 ). For Lorentzian triangulations, the 
analogous formula for a triangle with edge lengths i, t and / in a "traditional" coordinate 
system where i is the unit vector in the direction and / is the unit vector in the 1 
direction is: 



Note that this satisfies det(g) = —<iA 2 L (l,i,t). Again, we note that the metric in equa- 
tion (3.2) switches signature if the triangle inequalities are violated, however equation 
(3.3) remains in Lorentzian signature for all real values of i, t, I. In Euclidean Regge 
calculus the deficit angle at a vertex is given by 5 = 2n — J2 S @s, where the sum is over 
the interior angles 8 S of each triangle s incident on that vertex. In this case the deficit 
angle is always a real number. By considering equations (2.9)-(2.11) we see that for 
Lorentzian triangulations such as in figure 2 the angles summed around any vertex will 
be 2tc — irj. So now the deficit will be 5 = irj which is purely imaginary. Therefore 
curvature will be described by a pure boost, and this will be important in developing 
higher derivative 1+1 theories which we discuss in the conclusions. For any specific 
vertex the deficit angle is in general a complicated function of 12 edge lengths. 
However, in the limit in which all time-like edge lengths become null, (see figure 4 and 
section 4.1.1), the deficit angle simplifies to [19]: 




(3.2) 




(3.3) 




(3.4) 



This limit is discussed in more detail in section 4.1.1. 
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Figure 4. Traditional triangular lattice: Immediate neighborhood of a specific vertex in the 
limit where all the timelike edges become null. Only the spacelike lengths are labelled. 

3.3 Metric: symmetric 

In terms of the "symmetric" version of the triangulation one proceeds as follows: As 
in section 2, choose an initial coordinate system in which the vertices at the ends of 
the space-like edge of length S lie at (-S/2, —S/2) and (—5/2, 5/2) respectively. The 
third vertex is assumed to be in the forward lightcone of these two vertices at a point 
(t,x). For this calculation the key result we need is: 

- (t + S/2f + x 2 = ~\ + Tl + T 2 2 ) . (3.5) 

Now choose the time-like basis covector to be E ^ = rj^u \(T\ + T 2 )^ = (— (t + S/2), x), 
and the space-like basis covector to be E 1 ^ = (0,5/2). The (flat) metric inside the 
triangle is given by: 

g»v = Vap E\& u = -E\E\ + E\E\. (3.6) 
Calculating these out directly gives: 

Soo = -{t + 5/2) 2 + x 2 = -I Qs 3 + T 2 + T|) , (3.7) 
goi = gio = \xS= 1 -(T 2 -T?), (3.8) 
9ii = \S 2 . (3.9) 
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So the metric in this "symmetric" coordinate system is: 

Note that this satisfies det(g) = —Aj j (S,Ti,T 2 ). In general the deficit angle is quite 
complicated, however in the limit in which all time-like edge lengths become null, (see 
figure 5, and the discussion in section 4.1.1), the deficit angle simplifies to [19]: 

/ c-2 c-2\ 

S = i V = ihy(^y (3.11) 



s 2 




Figure 5. Regular triangular lattice: Immediate neighborhood of a specific vertex in the 
limit where all the timelike edges become null. Only the spacelike lengths are labelled. 



3.4 Action 

We can rewrite equations (3.1) in their simplicial form as: 

Z = [v[g] e iS[9] ->■ Z LFT = [ dn[{l 2 }} e tSL[{l}] , (3.12) 
J J{i} 

S = A f d 2 x^^S L = \Y,M{h}), (3.13) 

J s 

where {/} denotes the set of all edge lengths in the triangulation, and s is an index 
for each triangle in the triangulation. Furthermore {l s } is the set of edge lengths for 
the triangle s (we will denote these {l s ,t s ,i s }) and <i/z[{/}] is the measure on the set of 
edge lengths which we discuss later in the paper. It is instructive to compare equation 
(3.12) to the Euclidean path integral of QRC: 

Z QRC = [ dn{{l}) X [{1}) e- s ^™ , (3.14) 
J{l} 

S E = \ (3-15) 
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where is an "indicator" function on the set of edge lengths which is 1 when every 

edge length satisfies the triangle inequalities and otherwise. We can Wick rotate 
equation (3.12), t — >■ it, to Euclidean signature to obtain: 

Zlft -> Z LFTWlck = [ dfj,[{l} A ] e^K'W + / dfi[{l}d e- s *^ , 

J{i}a J{i}^ 

= Z QRC +[ dfi[{l}^]e- s ^ . (3.16) 

Here {/}a denotes the set of edge-length configurations in the Lorentzian theory which 
satisfy the Euclidean triangle inequalities, and {Z}^ denotes the set of configurations 
which explicitly violate them (possibly only one edge). [{/}$] is a pure imaginary 
number, and thus the Wick rotated Lorentzian path integral differs from the Euclidean 
path integral of QRC by 

Z LFT> wick ~ Z QRC = I (3.17) 

J {in 

The fact that Euclidean and Lorentzian quantum gravity are not related by a simple 
Wick rotation has been known for over a decade (see for example [1] and [2]). Equation 
(3.17) explicitly demonstrates that in this approach finding the relation between the 
two is at least as hard as fully solving the theory itself. 

Taking the Zlft of equation (3.12) as primary, we still need some method akin to 
Wick rotation so that the integral in equation (3.12) converges. The way we do this is 
to Wick rotate the lattice cosmological constant: A — » iX to give 

Zlft,e^ I dfx[{l}} e- SL[{l}] . (3.18) 
J{i] 

We note that this is analogous to the method used in CDT [1]; all results in the end 
will have to be analytically continued back in A. Note that with this prescription 
in Euclidean signature one still integrates over configurations satisfying the (vacuous) 
Lorentzian triangle inequalities weighted by the Lorentzian area functional. Analyzing 
equation (3.18) will be the focus of the next section. 

4 Analysis of the 1+1 Path Integral 

In QRC, the issue of which measure on edge lengths, to use has been the 

subject of some contention [20, 21]. In this paper we will follow reference [20] and will 
consider measures of the form 

d»[{l}\ = I[[V d (s)rlldll. (4.1) 

s ij 
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Here Vd(s) is the volume of the d-dimensional simplex s, (for the time being d = 2), 
and a is a parameter of the theory related to quantum gravity in the continuum by 

(d+l)[(l-u;)d-4] 
° = ^ , I 4 - 2 ) 

where a; is determined by the norm chosen in the superspace of metrics: 

\\5g\\ 2 = f d'x (g(x)r/ 2 G^[g{x)-M 8 9a p. (4.3) 



The case uj = is called the DeWitt measure, and in d = 4 corresponds to a = 0. 
However, in d = 2 we have o = —3/2. There is no unique preferred way to choose u, 
therefore in this paper we consider two specific measures: 

1. the simplest possible case — the uniform measure o = 0; 

2. the DeWitt measure a = —3/2. 

4.1 1+1 Lorentzian Fixed Triangulations with Uniform Measure 

Consider a triangulation of the torus which has N triangles and n edge lengths arranged 
as in figure 2. We wish to investigate equation (3.18) with the measure given in equation 
(4.1) for a = 0: 

/n 
JJ d/ 2 e -AEf^ ; (44) 

i=l 

where Aj is the area of the j th triangle in the triangulation. (Note that 2n = 3N = 
6N s N t , so n = 3N s N t .) Consider the Lorentzian area function for one triangle with 
space-like edge length I 2 = z, and time-like edge lengths t 2 = x and t 2 = y: 

A(z, x,y) = -^z 2 + x 2 + y 2 + 2zx + 2zy — 2xy. (4.5) 

This can be rewritten, by expressing (x, y, z) in standard spherical polar coordinates, 
as 

A(r, 6, (j)) = T - yjl + sin 2fl(cos + sin 0) - sin 2 9 sin 20 = T - f(6, 0), (4.6) 

which has zeros at r = and (r, 9, (f>) = (r, n/2, vr/4). The limits of integration for time- 
like edges are [0, oo), however the integral diverges if we allow the space-like edges to go 
to 0. This makes perfect sense, as the triangles remain well behaved when the time-like 
edges are 0, i.e. in the limit in which they become null. However, when the space-like 
edge goes to 0, the triangle becomes degenerate. Therefore we need to regulate the 
integration in such a way that the degenerate configurations are avoided. The simplest 
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way to do this would be to place a lower limit on the space-like edge lengths; however, 
for our purposes it is more convenient to remove a thin wedge from the single triangle 
configuration space by only allowing the azimuthal angle 6 to take values in [0 



' 2 



e 



We can make use of equation (4.6) by changing coordinates in the configuration space 
to n-dimensional hyperspherical coordinates, where r is now determined by the edge 
lengths of all the triangles, s, by r 2 = J2 S x 1 + vl + z l ■ Equation (4.4) becomes: 



oo 

A 



dVt n ^ / dr r"" 1 e — F(Q "- l)r . (4.7) 



Here F{Q n _i) is a complicated function of the n — 1 hyperspherical angles which has a 
zero only when all of the angles are chosen such that the area of every triangle is zero; 
these can be avoided provided we apply a regulation to the hyperspherical coordinates 
that is analogous to the one described above. For the uniform measure, the exact form 
of F(Q n _i) is unimportant in determining scaling relations since we can write: 

r(n) (af(Fvo)" = a(£) ^W^' (4 ' 8) 

where we have used the topological fact 2n = 3N (each of a triangle's 3 edges is shared 
by 2 triangles). We have introduced a(e) which is just a numerical factor depending 
on the UV cutoff e given by: 



a ^ = J s / n ^mz^' (4 ' 9) 

where the zero of F(Q n _i) is avoided by only integrating over the e regulated surface 
S e . The numerical constant is unimportant in calculating how the average area of a 
triangle scales with A: 

(A) _ 1 d\n(Z) _ 3 1 

N N dX 2 X { ' ' 

This scaling relation for the area is identical to that found in QRC [22], this is essentially 
because in both path integrals, every edge length can be rescaled by I 2 —> l 2 /X. Thus, 
for the uniform measure, one must look beyond the average area scaling relation to find 
differences between the two theories. 

Note also that in view of the fixed topology, (of the triangulated manifold), and 
the Gauss-Bonnet theorem, we have 

J2 S v = 0- (4-11) 

V 

This is a sum over all vertices v, and is guaranteed to vanish at the kinematical level 
for each individual geometry in the path integral. (It is very easy to verify this result 
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explicitly for the limit where all timelike edges become null.) This implies an exact 
result for the spacetime average value of the Ricci scalar which holds for any choice of 
the measure: 

(R) = 0. (4.12) 
4.1.1 Causal Diamond Truncated Model 

We now discuss a simplifying truncation of the model which will also be of use in the 
next section, but is worth considering here: that is the limit in which all time-like edges 
become null. In this case the triangulation is made up of "causal diamonds" formed by 
the two triangles connected by each space-like slice; these causal diamonds have area: 

A(s) = \s\ (4.13) 

With reference to equation (3.10), the metric inside each causal diamond becomes 

s 2 f-1 0\ s 2 



^ = T V o i) = J ^' (414) 

where s is acting like a conformal factor. With this truncation the integration in 
equation (4.4) is trivial: 

// roc \ N/2 yN/2 

n<*(«?)e-* AE ^=(jf dxe-^j =jm . (4.15) 

(There are N/2 = N s N t spacelike edges.) This yields the simplified scaling relation 

It is also trivial to find the average length of the spatial slices in this truncation: 
(L) = N s (s)=N s (^J^dxx 1/2 e-^^dxe-^ = N s ^ X' 1/2 . (4.17) 

Note that is is now very easy to generate spatially large universes — one just has 
to choose N s (and A) appropriately. Geometrically each of these null-edged triangles 
extends a distance s/2 into the time direction, so the spacetime has temporal extent 

<^*f^^ < 418 > 

So it is also very easy to generate long-lived universes — one just has to choose Nt 
(and A) appropriately. Now observe 

(A) (s 2 ) 4 



<L>m w -- 1 ' 27 - (419) 

This can be interpreted as saying that the "typical" geometry is not too crumpled, and 
suggests we are dealing with a relatively "smooth" phase. 
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4.2 1+1 Lorentzian Fixed Triangulations with the DeWitt Measure 

The analysis in this section proceeds in an identical way to the previous section, except 
now with a non-trivial measure, specifically a = — |. The path integral is now: 

/N n 
n^]~ 3/2 IW ' A>: ' :1 '. (4.20) 
j=l i=l 

We can again switch to hyperspherical coordinates, where now we denote (9j,4>j) as 
the hyperspherical angles for the edge lengths corresponding to the triangle j. With 
this we can write: 

Z = jtl [iM^i)]' 372 ^ 1 e-^Uf^) drdtln^ 



3=1 

j AT 



dO n _i 



G(Q n _!)3/2 



r -l+n-3N/2 e -±F { n n _ l)r ^ (4 21) 



where G(f2 n _i) = FT/ /(fy> $7) anc ^ ^(^n-i) = fi^jt^j)- Both these functions 
are more complicated than just products and sums of functions of the form f(6,4>) 
in equation (4.6), due to the fact that the same angles will appear multiple times 
depending on the incidence matrix of the chosen triangulation. We also note that 
F(O n _i) < y/2N and G(fi n _i) < 2 7V/2 . We again make use of the relation 2n = 3N to 
write: 



Z = 8 N I / r- 1 e-^nn-i)' dr 



G(n n _!) 3 / 2 

G(fi n _!) 3 / 2 



oo 

A 



0; 



(4.22) 



In contrast to equation (4.7) to ensure convergence we have now had to introduce 
an explicit ultraviolet cutoff in the r integration, as well as the angular regulation 
introduced in section 4.1. This means that the smallest the space-like edge can become 
is esin(e) ~ e 2 . Here r[0;x] is the incomplete Gamma function which has the series 
expansion: 

(-l) k+1 x k 



r[0;x] = -IE - In(x) + l kl • ( 4 - 23 ) 



k=l 



This converges for all x e (0, oo). Equation (4.20) is then given by: 



Z = a(e) + /3(e) In (^j + a k (e) e k X k , (4.24) 
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with 



a(e) 
a k (e) 



N r gu r / 4 
] St \ 



- IE 



dtt 



n-l 



,8. G(tVl) 3 / 2 ' 

Js c 



A k kk\ 



n-l 



G(a„_o 3/2 ' 



(4.25) 
(4.26) 
(4.27) 



Here a(e), /3(e), and ajt(e) are numerical factors that depend on the ultraviolet cutoff 



e. 



In contrast to the situation for uniform measure, the scaling relation for the average 
area is now non-trivial: 



(A) 



1 



X V«(e) - /3(e) ln(A) + ZZi «fc(e)e fc A^ 



(4.28) 



While this relation demonstrates that the case of the DeWitt measure is less trivial 
than the case of the uniform measure, it does not give much more information than 
that. It may be possible that for a fixed value of e there is a value of A where the 
average area vanishes or diverges, however (without knowing the precise values of the 
cutoff-dependent quantities a, (3, and a*.) determining if this can happen is impossible. 
Attempts at performing the integrals in equations (4.25)-(4.27) will run into similar ob- 
stacles to those presented in the next section, therefore ultimately a numerical analysis 
will be necessary for further investigation of the full model. 

4.2.1 Causal Diamond Truncated Model 

To further investigate equation (4.20) we can look at it in the causal diamond truncation 
outlined in section 4.1.1. Again, the integration again becomes tractable, but still 
exhibits behavior similar to equation (4.22): 

d(sf) e"^*^ 




(4.29) 
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(A) 



Figure 6. A plot of -jt 
DeWitt measure and e 



lis 2 ) and & 



(s) for the causal diamond truncated model with 



0.1. 



The average area then scales as: 



{A) - 1 r &'l - 1 ./^ + 0(f) . (4.30) 



iV 4A AT" „_ Af / 2 nn.Aei 4A 



2_ p-Ae/2 _ p [1. Ae] 
\e L2' 2 J 

This is a limiting case of equation (4.28) and is well behaved for all finite A (see figure 
6), which suggests that equation (4.28) is as well. The average length of the spatial 
slices can also be computed by using (L) = N s (s), while the average temporal extent 
of the spacetime is (T) = ^N t (s), where it is easy to establish that 

W - . ..^1' - 4 r- (D - ^ + Oi^ny (431) 



2e-¥-V2iAr[i;f] 2 V 
(See figure 6.) In this situation we furthermore have 

(A) (O 2r[i;f]( V ^e-^-r[|;f] 



<L> (T) ( a >» r[0 



■ eAl2 
' 2 J 



(4.32) 



(See figure 7.) As Ae -> this blows up as O ((Ae) -1 / 2 ^ + ln(Ae/2)]" 2 ). In contrast 
as Ae — > oo this smoothly tends to unity from above. 

{A) { ' ,2} 1 + (4-33) 



(L) (T) (sy- (Xe)' \(Xe) 

There is no sign of any phase transition, and as Ae — > oo becomes arbitrarily smooth. 
Again it is very easy to get arbitrarily large universes, both in spatial and temporal 
extent. 
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X 
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Figure 7. A plot of tA -pj* = -he? for the causal diamond truncated model with DeWitt 
measure and e = 0.1. 



4.3 Summary 

For both the uniform and DeWitt measure, the path integral seems to be well behaved 
in A, meaning that it does not exhibit critical behavior. (There do not seem to be any 
phase transitions.) The non existence of a critical point means that this model does 
not seem to have the usual continuum limit based on a 2 nd order phase transition. The 
non-criticality of the pure 2D gravity QRC path integral was shown in [14], however 
once a higher derivative R 2 term was added to the action and the partition function 
was restricted to fixed area critical behavior was observed numerically [15]. Future 
numerical work might examine the analogous behavior for this Lorentzian model with 
an R 2 term in the action (see [23, 24] for implementations in QRC and DT). As we 
mention in the discussion, there are reasons to believe the absence of a phase transition 
might actually be beneficial. 

5 Loop-to-Loop Amplitudes 

5.1 Ai-step Amplitude 

An important quantitiy we would like to be able to calculate in the model is the loop-to- 
loop amplitude, G\(L, K; T), which is the amplitude for a loop of length L to propagate 
to a loop of length K in "time" T. It will be a sum over all interpolating geometries of 
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K 




Figure 8. A Loop-to-Loop Geometry 



the type shown in figure 8. This amplitude has been calculated for CDTs in [1] and in 
a continuum calculation in proper-time gauge in [25]. 

Unlike CDTs, in the current model there is no natural way to define the T appearing 
in the loop-to-loop amplitude G\(L, K; T). There is however a natural way to use the 
fixed triangulation to define the "iV r step" loop-to-loop amplitude: 



This amplitude will be given by equation (3.18) for a triangulated cylinder, where 
the spatial boundary lengths are fixed to be L and K. For an N t x N s triangulation 
of the cylinder in the form of figure 2, let the space-like slices have length Li with 
i — 0, 1, N s . The amplitude (5.1) (for the uniform measure) is then given by: 



G\(L, K; N t ) = / Y\ dL? G\(L, L\\ 1) G\(Lx, L 2 ; 1) ... G x {L Nt ^,K;l), (5.2) 



where G\(Li, Li+i; 1) is acting as a single-step temporal transfer matrix which we will 
call the strip amplitude. To construct the strip amplitude, consider the triangulated 
strip shown in figure 9. The space-like edges are fixed to be = Li/N s and = 
Li + i/N s , however the vertical time-like edges labelled i n , and the diagonal time-like 
edges labelled t n with n = 1,...,N S , are dynamical degrees of freedom. 
The strip amplitude is then given by: 



G\(L,K;N t ). 



(5.1) 




8=1 



G x (L h L i+1 ;l) = / \[dtlde n e 




■N s 

'71 — 1 



A l (li ,tn itn — 1 ) + (k+1 jtn fin ) 



(5.3) 
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Figure 9. A Triangulated Strip 



This can now be rewritten in terms of single-step spatial transfer matrices, U and V, 
by making the identification x 2n -i = t^, x 2n = t^, under which equation (5.3) becomes 

POO 

G\(Li, L i+ i, 1) = / dxidx 2 . . .dx 2 N s 
Jo 

xU(x u x 2 ) V(x 2 ,x 3 ) . . .U{x 2Ns -x,x 2Ns ) V{x 2Ns ,xi) 
= Tr[(UV) Ns ], (5.4) 



where 



V (x n ^ x r< 



—^A L {li +1 ,x n ,x m ) 

j 

J -XA L (k,X n ,X m ) 



(5.5) 
(5.6) 



Thus, finding the strip amplitude is reduced to finding the spectrum (in fact since we 
can take N s to be large, the largest eigenvalue) of the integral operator 



UV(x, y) = I dz e 
'o 



-f (v / ^+ a ) 2 + fo2 +v / ( 2 + c ) 2+d2 ) 



Here 



(5.7) 



a — x; 
b = 2l i+ iy/x; 



d = 2l iy /y. 



(5i 



It is unclear if any explicit solution to equation (5.7) can be found in terms of special 
functions. However in the special case where U = lj, the diagonal of UV can be readily 
computed: 



UV(x,x) 



b I < - 

'sinh -1 (a/6) 



2ZiVx Ki 



^ cosh W cosh(t) dt 
If-x 



\li\fx , sinh 



-1 / "i 



(5.9) 



Here Ki[x,y] is an incomplete modified Bessel function of the second kind [26, 27]. 
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Despite the difficulty of the integral appearing in equation (5.7), we can make 
progress by investigating the integral operators U and V with kernels defined in equa- 
tions (5.5) and (5.6) individually. Considered as operators on Lj|(0, oo), that is the 
space of all real square integrable functions on [0, oo), they are positive, symmetric and 
bounded. For purposes of the following calculation temporarily identify I 2 , l 2 +1 — > I, 
and A/4 — >■ A, and use the inequality 



+ x + y) 2 - Axy < I + x + y, (5.10) 

which implies 

e -\y/(l+x+y) 2 -4xy > e ~X{x+y+l) ^ (5-11) 

It is now easy to show that U and V are positive definite with respect to the L 2 inner 
product: 

/oo /roc 
dxf(x)lj dyU(x,y)f(y) 

oo poo 



dxdy e- x ^ {l+x+y)2 - Axy f{x)f{y) 



JO 

oo poo 



> / / dxdye- x « +x+ rtf{x)f(y) 
Jo Jo 



XI I I J„ „— Xx . 



dxe~ Xx f(x)j >0. (5.12) 

In addition the trace of U and V are finite, so the operators are trace class, and the 
traces are given by: 

POC poo cyf-i _|_ \j2\ xl 2 

Tr[U] = / dxU(x,x) = / dxe- x ^+ ilx = 1 ^ J e"^, (5.13) 

Jo Jo A ij 

and 

Tr\V] = 1 e— (5.14) 

These properties of U and allow us to apply the Hilbert space generalization of the 
trace inequality [28, 29]: 

Tr[(AB) M ] < (Tr[A]) M (Tr[B]) M . (5.15) 

Thus, the strip amplitude is bounded above by: 

G x {L ii L i+1 ;l)=Tr[{UV) N '] 

< (Tr[U}) N ° (Tr[V}) Ns 

N (1 + Xlf) Ns (1 + Xlf +1 ) Ns AiV s , ;2+/2 . 

= 4^ ^ ^ — ^ * +lJ e—^b+W. 5.16 

Note that equation (5.16) is well behaved for all finite values of A. 



- 20 - 



5.2 Proper Time Amplitude 

As we mentioned above, there is no natural "proper time" parameter appearing in the 
model. We can however introduce one, and the most obvious way to do so is to fix 
the "vertical" time-like edge lengths of the triangulation of figure 2 to a fixed common 
value r = T/N t . This can be done by fixing all t n = r in figure 9, while allowing 
the t n to remain dynamical. Once this is done, the flat space (Minkowski) solution 
for an N t x iV s triangulated cylinder with time-like length (height) T and space-like 
length (radius) L is clearly -t 2 = -r 2 + I 2 = -(T/N t ) 2 + (L/N s ) 2 , where r = T/N t is 
the time-like edge length of the vertical edges, I = L/N s is the space-like edge length 
of the horizontal edges, and t is the time-like length of diagonal edges. Non-flat 1+1 
geometries correspond to allowing t to fluctuate from its Minkowski value. As above, 
to calculate the loop amplitude G\(L, K; T) one must integrate over the diagonal edge 
lengths: 

G X (L,K;T)= / Y[ dLlGl{L,L X ]l) Gl(Lt,L 2 ;l) ... Gl(L Nt ^, K;l), (5.17) 
i=i 

where now we need to find the strip amplitude with fixed "vertical" time-like edges 
G T x (Li, L i+ i, 1). Notice that when the vertical edge lengths are fixed, the dynamical 
edge lengths only couple two neighboring triangles together, thus the strip amplitude 
factorizes: 

(/•oo \ N 3 

J dxe- XA WW + ^ 2 n , (5.18) 

where 
A(x; a, b,c) — — 

(V^ 2 + 2(a - cjx + (a + c) 2 + x 2 + 2(6 - cjx + {b + c) 2 ) . (5.19) 

Again, the integral found in equation (5.18) is of the same form as that in equation 
(5.7) and whether a closed form solution can be found is at this stage unclear. However, 
the diagonal elements of the strip amplitude can be found by solving the integral in 
equation (5.18) with Lj = L i+X = L: 



2LT „ T ALT . . ! f(N s /N t )L 2 - (N t /N S )T 2X 1 ^ V 



iVVVV Sinh ~ V 2LT 



(5.20) 



The strip amplitude can also be estimated for large A by performing a saddle point 
approximation. In this case, with N t = N s for simplicity, the strip amplitude is ap- 
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proximated by: 




3 \ AT s /2 



x e 4N ' 



A 



^^/(Li-L i+1 ) 2 +4T 2 (Li+L i+1 ) T £ rp2. 



G\(Li, L i+ i; 1 



) = 



(5.21) 




v 



5.3 Summary 

The expressions derived above can be compared with the analgous expression in CDT 
[1]. The most important fact is that the CDT expression has a non-zero critical value 
for A, while we have shown that equation (5.4) is bounded above by an expression 
which is well behaved for all finite A. It should not be unexpected that the quantities 
calculated in the different theories should be different, perhaps radically different, at 
the discrete level. However, it is clear that while in CDT a continuum limit can be 
taken by looking at a scaling limit near the critical value of A, this does not seem to 
be possible for the current Lorentzian model. The only sensible way to compare the 
two theories is in their continuum limits, thus it is important that some way of taking 
a continuum limit be found for this model. As discussed in section 4 a possible way to 
do this would be to introduce an R 2 term into the action. That will be the topic of 
further research. 

6 The non-occurrence of spikes in LFT 

An problematic issue with Euclidean QRC, raised in reference [21], is what can be called 
the "occurrence of spikes" — wherein "spiky" geometries dominate the path integral. 
Specifically, in (Euclidean) QRC for a fixed value of space-time area A there exists a 
finite n such that for any edge length in the triangulation 



We will now argue that in the full LFT model these spikes do not appear, and we 
will show by explicit calculation that they also do not appear in the Causal Diamond 
truncation. 



(l n )A = OO. 



(6.1) 
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Figure 10. A hexagonal neighborhood of a vertex in the triangulation. In the Euclidean 
case spikes appear at the central vertex when all of the incident edge lengths become large 
L, while all of the boundary edge lengths become small, I = j in order to keep the total area 
finite. This situation cannot occur in the Lorentzian case. 

Because the triangulation introduced for our model is hexagonal (see figure 10), 
we will focus on the analogous geometry in reference [21]. In the Euclidean model it is 
demonstrated that spikes appear by considering regions of the configuration space in 
which all of the incident edge lengths on some vertex become very large ~ L. To keep 
the total area of the triangulation fixed all of the edge lengths of the boundary hexagon 
become correspondingly very small ~ / = 1/L. Due to the form of the Euclidean 
area function (2.7), one can take the limit L — > oo and keep the total area of the 
triangulation fixed since the area of each triangle in the hexagon scales like 



However, this cannot occur in the Lorentzian model because due to the form of the 
Lorentzian area function (2.6) if one tries to take the same limit the now positive quartic 
terms dominate in four of the hexagon's six triangles: 



In two of the hexagon's triangles (the ones with both time-like edges ~ L) the quartic 
terms can be arranged to cancel and the areas remain bounded, however the total area 
of the hexagon must diverge. In fact, just one of the triangle areas diverging is enough 
to show that the hexagonal spikes that occur in Euclidean QRC cannot occur the LFT 
model. 

In fact, due to the form of (2.6) the space-like edge lengths S of the triangulation 
can never go beyond a finite bound if the area is to remain fixed. As alluded to above, 
the only way one can scale the edge lengths of a Lorentzian triangle such that some 
edge lengths diverge while the total area remains fixed is to set 7\ = T 2 = L and 



A E ~ IL = 1. 



(6.2) 



A L ~ L 2 ->■ oo. 



(6.3) 
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S — 1/L, in which case the area scales like 

A L ~ 1. (6.4) 

Note however that for the triangulation we have considered in this model, if we attempt 
to scale a triangle in this way then the areas of the adjacent triangles, that is, the ones 
to the left and right in the strip (see figure 9), will diverge unless they are scaled in 
an identical manner. This argument then extends to every triangle in the strip which 
contains the original triangle. Thus these regions in configuration space correspond to 
"pinched" geometries, that is geometries in which the spatial extent of a strip shrinks 
to zero while the temporal extent diverges (see figure 11). 



Figure 11. A Pinched Geometry. Compare with the "A" phase of CDT. 

The extent to which these pinched geometries will contribute to the partition func- 
tion will depend on the how the parameters of the theory are tuned. We note that 
geometries very similar to this are dominant in the so-called "A" phase of CDT [30]. 
It is possible then, perhaps only after adding an R 2 term to the action, that for fixed 
space-time area the general LFT model has phases where the pinched geometries are 
dominant and others where they are not. 

At any rate, these pinched geometries can be automatically regulated away by 
going to one of the two truncations introduced in the previous sections. In the proper 
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time truncation, the pinched geometries clearly cannot appear since half of the time- 
like edges in the triangulation are by definition frozen to a fixed value. This means that 
the scaling procedure outlined above cannot be performed. Perhaps more interesting is 
the Causal Diamond truncation in which the quantity (l n )A can be calculated exactly. 

In the Causal Diamond truncation with fixed space-time area we need to consider 
the partition function 



Z[A] = I Y[dXi etXi* dt^Yl 

i \ i 



cyM 



(6.5) 

where M = N/2 = N s N t is the number of space-like edge lengths and Vm-i is the 
volume of a right-angled M — 1 simplex (in configuration space) which has all of its 
right-angle edge lengths equal to 2A. Then 

Z(r) A = 2e~ XA dxi--- dxy I dx M -i{x M -i)^ 

Jo Jo Jo 



— — / dxi ■ ■ ■ / dx M -2 I 2A - ^ Xj - 

■ + 1 JO Jo \ „■_•, 



= 2 e XA n ^ 2 - 1 J dxi-- - I dx M -2 \2A- > \ x { - x M -2 

2 2M-l+f ^M-l+f e -\A 

~ (n + 2)(n + 4) • • • (n + 2(M - 1)) ' 

One performs the iterated integration above by repeatedly making the linear change of 
variables x = 2A — x i ~ x j- Thus the expectation value is given by: 

{l n ) A = 2^ 1 n T{M)T ^} At. (6.6) 

w r(M + f) v ; 

This is finite for any finite n and thus the appearance of spikes is completely removed 
from the Causal Diamond model. A quick sanity check with n = 2 reveals: 

exactly as should be expected for a fixed area. Note also that 

(l n ) A ~ A? ~ (Z>S, (6.8) 
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so a variant of the scaling argument presented in the discussion following equation (37) 
of reference [21] implies that the fractal dimension for the Causal Diamond model is 
simply its topological dimension, namely 2. We remark that in the "thermodynamic 
limit" , where M — > oo holding n and A fixed, we have the stronger result 

(Ha = r(l + |) (l 2 )T , (6-9) 

while the average area per triangle in this limit is simply \{l 2 ) a- 

In view of these comments, it is therefore clear that (at least in two particularly 
interesting truncations of the LFT model) the spikes which plague the Euclidean QRC 
simply do not occur. We have also argued that the QRC spikes cannot occur in the full 
model either, however there are still regions of the configuration space which correspond 
to "pinched" geometries. Whether these geometries will dominate the partition function 
(and if they do whether or not that would be pathological or interesting) will be the 
topic of further research. 




Figure 12. (3,1) tetrahedron Figure 13. (2,2) tetrahedron 



7 Higher Dimensions 

The model discussed in the previous sections can be extended into higher dimensions 
with a few modifications. The first is that in n dimensions we now need to consider 
the Einstein-Hilbert term in the gravitational action [31]: 

S =\f d n x ^{R-2K) ^ KY,V h 5 h - xY^Vl^s), (7.1) 

^ h s 

where h denotes the n — 2 dimensional "hinges" in the triangulation along which cur- 
vature is concentrated, Vh is the Lorentzian volume of a hinge, 5h is the deficit angle of 
the hinge and Vl(s) is the Lorentzian volume of the n-dimensional simplices s. 
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This model will differ from QRC in three ways: first, the volumes of the simplices 
will be determined by their Lorentzian formulas, which can be derived from the anal- 
ogous Euclidean formulas in a way similar to equation (2.8); the exact procedure is 
outlined in the next section. Second, the nature of the deficit angle of a hinge will 
depend on whether that hinge is space-like, or time-like, or perhaps even null (see [13]). 
Finally, there will be generalized triangle inequality constraints on the edge length 
configuration space, however these will not be the constraints found in the Euclidean 
theory as there will be certain simplices in the triangulation whose edge-lengths are 
constrained and others whose edge-lengths are not; we will derive the Lorentzian gen- 
eralized triangle inequalities in 2 + 1 dimensions below. The Euclidean and Lorentzian 
QRC models will be related in a manner similar to equation (3.16), which is clearly 
much more involved then a simple, naive Wick rotation. In light of this, studying 
the Lorentzian QRC calculus numerically will involve adapting previous studies of Eu- 
clidean QRC by implementing the changes described above. Analytic investigations 
will be mainly constrained by the same obstacles encountered in Euclidean QRC, and 
will therefore by strongly limited. 

In developing the model in 2+1 dimensions it will, in general, be impossible to 
triangulate the manifold in such a way that every face of the Lorentzian tetrahedron 
can be taken to be a Lorentzian triangle. There are 4 types of Lorentzian tetrahedron; 
two of them are very simply related to the Euclidean tetrahedron: one with all six edges 
space-like, i.e., the standard Euclidean tetrahedron itself with the usual Euclidean 
triangle inequalities, and the trivial generalization in which all six edge lengths are 
time-like, in this case all triangle inequalities are reversed (due to the —1 in the metric) 
but the tetrahedral conditions are unchanged. 

As described in [32], for the types of manifolds we are interested in, there are two 
types of tetrahedra which need to be considered: the (3,1) type in which a Euclidean 
triangle in one time slice is connected to a vertex in a subsequent time slice by three 
time-like edges (see figure 12); the other is the (2,2) type in which there are two space- 
like edges, each in a separate time slice, connected by four time-like edges (see figure 
13). As we will soon demonstrate the (3,1) tetrahedron is analogous to the Lorentzian 
triangle in that its time-like edge lengths are completely unconstrained (its space-like 
edge lengths form a Euclidean triangle and therefore must satisfy the Euclidean triangle 
inequalities); however, all edge lengths of the (2,2) tetrahedron must satisfy complicated 
constraints. 

7.1 (3,1) Lorentzian Tetrahedron 

The question is, under what conditions is the sextuplet S^i) = {Si, £2, S3, T±, T 2 , T 3 }, 
where Si is opposite T i} a (3,1) tetrahedron? First the tetrahedron has four faces, one 
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is a standard Euclidean triangle with edge lengths Si, S 2 and S 3 , and the remaining 
three are Lorentzian triangles with one space-like edge and two time-like edges. Thus 
the only requirement on the faces is that the edges of the space-like triangle satisfy the 
Euclidean triangle inequalities: 



51 < S2 + S3, 

52 < S3 + Si, 

53 < Si + S2. 



(7.2) 



Note that having the faces satisfy these conditions only means that the 4 triangles 
formed by S(s,x) are at this stage merely a "facial net" in the language of [33]. To form 
a Lorentzian tetrahedron the volume must be a positive real number, equivalently the 
Cayley-Menger determinant must be positive. 

The volume of a Euclidean n-simplex is given by: 



'I 



"I 



\n+l 



2 n (n\ 



(7.3) 



where Efj is the (n + 2) x (n + 2) matrix of squared Euclidean distances between the 



n + l vertices i and j in the simplex, augmented with an additional row and column 
defined by E^ = 0, E 2 = 1 and E^ = 1. The determinant of this matrix is known as 
the Cayley-Menger determinant. For example, a triangle with side lengths a, b, and c, 
has area given by: 



1 

16 



111 

1 a 2 b 2 
1 a 2 c 2 
1 b 2 c 2 



(7.4) 



which results in Heron's formula (2.7). 

Note that this determinant is positive iff (if and only if) a, b and c satisfy the 
triangle inequalities. To find the analogous formula for Lorentzian simplices we make 
use of the following observation. Given a Euclidean simplex with edge lengths {Li}, if 
we partition this set into two sets {Si} and {T{} then the analogous Lorentzian simplex, 
with time-like edge lengths {Tj} and space-like edge lengths {S,}, has volume given by: 



v L [{T i } 7 {s i }] = -iv E [{m},m}. 



(7.5) 



Thus, Vl 



— V|, and therefore the volume of a Lorentzian n-simplex is given by: 



(-1)" |f2 | 
2 n( n f)2l ij\ 



(7.6) 



- 28 - 



where L 2 - is the now Lorentzian distance squared from vertex % to j. Thus the area of 
a Lorentzian triangle with space-like edge length S and time-like edge lengths T\ and 
T 2 is given by 



(A 2 )l 



1 

16 



1 


S 2 



1 1 

2 rp2 



S 2 —T{ 


-T| 



-T 2 

1 2 



(7.7) 



which results in the formula we derived before, equation (2.6). Note that this determi- 
nant is positive for any choice of S, T 1; T 2 . 

Finally, we come to the volume formula for the (3, 1) tetrahedron with edge lengths 
as above: 

1 


s 2 



i) 
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si 
si 
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-T| 
-Tl 
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Remarkably (7.8) is positive for arbitrary real values of T 1; T 2 , and T 3 , provided that 
(7.2) is satisfied. That is, one need only apply triangle inequalities to the single Eu- 
clidean face, the three Lorentzian faces are unconstrained. A symmetric way of writing 
out the volume is: 



Vol(3, 1) 



1 

12 



S\T 2 X (S 2 + S 2 - Sf + T 2 - T 2 - T 2 ) 



+ S 2 2 T 2 2 {Sl + Sl - S 2 2 + Tl - Tl 



+ S 2 3 Tl (Sl + S 



si 



rp2 

1 3 



Tl 



Tl) 
Tl) 



+ - 2 (SI 



Tl){S 2 2 + T 2 2 ){Sl + Tl) + -(Sl - Tl){S 2 2 - Tl){S, 



Tl) 



1/2 



(7.9) 



Another way to write out the volume, with x — ^(S 2 + S 3 — Si), y = ^(S 3 + Si — S 2 ) , 
z = l(Si + S 2 - S3), and s = ^(Si + S 2 + S3) = x + y + z, is: 



Vol(3, 1) 



1 



12 



a - x) 2 Tl + (s- y) 2 n + (s - z) 2 T^ 

+ 2{s - x) 2 Tl(sx - yz) + 2{s - y) 2 Tl(sy - xz) + 2(s - z) 2 Tl(sz - xy) 
- 2TlTl(sz - xy) - 2T 1 2 T 3 3 ( S y - xz) - 2T 2 2 T 3 2 ( S x - yz) 



+ {s-x) 2 {s-y) 2 {s~zf 



1/2 



(7.10) 
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The formula for the cosine of the dihedral angles is found by making the appropriate 
substitutions into the formula for a Euclidean tetrahedron. For the Euclidean case, with 
edge lengths a opposite to d, b opposite to e, and c opposite to /, and the four faces 
given by Aabf, Abed, Aace, and Adef, the cosine of the dihedral angle at the edge a 
is given by: 



COs6> a 



[-a 4 



e 2 )(b 2 - f 2 ) + a 2 (b 2 + c 2 - 2d 2 + e 2 + f 2 )} 



(7.11) 



16A E (a, c, e) A E (a,b,f) 
In this case then the cosine of the dihedral angle at a space-like edge cos 9g and the 
dihedral angle at a time-like edge cos 9t are given by: 



Sl ~ 16iA L {S 1} T 3} T 3 )A E {S 1} S 2} S 3 ) ' U } 

[Tf + (S 2 + T 2 )(S 2 + Tj) - T 2 (T[ + T 2 + 2S 2 - S 2 - S 2 )] 
C ° S9Ti ~ 16 A L (S 2 ,T U T 3 ) A L (S 3 ,T1,T2) " (7 ' 13) 

The sine of the dihedral angle is much easier to find by using the Lorentzian analog 
of the Euclidean <i-dimensional sine law: 



d 



V n V r 



where V n is the volume of the n-simplex, V^_2 is the volume of the hinge at which the 
dihedral angle is being considered, and V n -\ and V^_ 1 are the volumes of the two faces 
which share the hinge. From this we can derive the sine of the dihedral angle about a 
space-like edge length, sin^s, and about a time-like edge length, sin^-r: 

35*1 V(3 ) i)(S'i, S 2 , S , 3 ,Ti,T 2 ,T 3 ) 



sin 9 



2 A L (S 1 ,T 2 ,T 3 )A E (S 1 ,S 2 ,S 3 ) 
3Ti V( 3i i)(S'i, ^2, S , 3 ,Ti,T 2 ,T 3 ) 



(7.15) 

2 A L (S 2 ,T x ,Ts) A L {S 3 ,T u T 2 y (7 ' 16) 
We note that equations (7.10), (7.12), (7.13), (7.15), and (7.16), reproduce the results 
found in [32] for the special case S\ = S 2 = S 3 = 1 and T\ = T 2 = T 3 = y/a. 



Si 



sin 9t x 



7.2 (2,2) Lorentzian Tetrahedron 

The sextuple of edge lengths for the (2,2) tetrahedron is 5(2,2) = {Si, S 2 , T\,T 2 , T 3 , T4}, 
where Si is opposite S 2 , T\ is opposite T 3 and T 2 is opposite T 4 . The volume is then 
given by: 

1 1 1 1 

1 Si -T 2 
1 S 2 



3^(2,2) 



1 

'288 



_rp 2 

-n -if 



1-T 2 -T 2 S 2 



-T 2 



-Tf 



5.? 







(7.17) 
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The constraints on the edge lengths are now much more complicated than simple tri- 
angle inequalities. If the edge lengths are considered as vectors, then one can derive 
T\ + (— T 2 ) + T 3 + (— T4) = 0. That is, starting at the {S\,T\,T±) vertex in figure 
13 one can follow the closed path defined by the previous equation. This implies 
that Ti + T 3 = T 2 + T 4 , which is equivalent to saying the four time-like edges form a 
non-planar quadrilateral, and are therefore constrained — they cannot take on arbi- 
trary values. 1 An example of a sextuple which is an impossible (2,2) tetrahedron is 
£(2,2) — {1/2, 1, \/2, 1/V2) 1) 1}- A symmetrical representation of the volume formula 
is: 



Vol(2,2) = 1 



(S 2 + S 2 + T 2 + Ti + T 3 2 + T 2 ) 

-T 2 T 3 2 (S 2 + S 2 2 + T 2 - T| + T 3 2 - T 2 ) 
-T 2 T 2 (S 2 + S 2 2 - T 2 + T| - T 2 + T 2 ) 



+ \{Sl + Sl){Tl + Tl){n + Tl) 
-\{Sl-Sl){Tl-Tl){Ti-Tl) 



1/2 



(7.18) 



It appears that within the framework outlined in this paper the (2,2) tetrahedron re- 
introduces the difficulty of configuration space constraints on the time-like edge lengths. 
This will impede analytic calculations in the same way the Euclidean constraints impede 
calculations in QRC Therefore studying this model in higher dimensions will best be 
done numerically. 

In the same manner as in section 7.1 we can find the cosine of the dihedral angles 
at a space-like edge length, cos 9s, and a time-like edge length, cos 9t- 

[Sf - (T 2 - T|)(T 2 - Tj) + Sf (T 2 + T| + Tf + Tj + 2S 2 )] 
Sl ~~ IQA L {S U T U T 2 )A L {S X ^T A ) ' {LU} 

n P? ~ (Sf + Tj){Sl + T 2 ) + T 2 (Sf + SI + 2T| - T 2 - T 2 )] 
COS ^- 16A L (SuTuT,) A L (S 2 , Tl ,T 4 ) " (7 ' 20) 

We can also find the sine of the dihedral angles at a space-like edge length, sin 9g, and 
a time-like edge length, sin 9t'- 

. Q .33i V(2,2){Si,S 2 ,Ti,T 2 ,T 3 ,T 4t ) 

" -*-2T AUSuTuTJ A L (S U T 3 ,TA> (? ' 21) 
_ 3Ti V( 2j2 )(S'i, S 2 ,T 1 ,T 2 ,T 3 ,T 4 ) 

sm ^ Tl " ~2~A L (,Si,Ti,r2) A L (s 2 ^ Ti y (7 ' 22) 



1 Though we can easily solve this particular constraint by setting T1.3 = a ± 6 and T2.4 = a ± c, this 
merely moves the problem elsewhere. We then have to constrain a, b, c by demanding that a ± 6 and 
a ± c all be timelike. 
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Equations (7.18)-(7.22) also reduce to the values found in reference [32] for the specific 
values S\ — S2 — S3 — 1 and T\ = T 2 = T 3 = ^/a. 

7.3 3 + 1 and Beyond 

The procedure carried out in sections 7.1 and 7.2 can be applied in exactly the same 
way in any dimension n — d + 1, however with increasing difficulty as the number of 
edge lengths grows like n 2 . In 3 + 1 dimensions there are again two types of 4-simplex 
that need to be considered: the (4,1) type which is the analog of the (3,1) tetrahedron, 
and the (3,2) type which is the analog of the (2,2) tetrahedron [32]. Finding the volume 
and the sine of the dihedral angles is simply a matter of applying equations (7.6) and 
(7.11) respectively. Finding the cosine of the dihedral angle is more involved but can 
be done using the higher dimensional analog of the cosine law. Explicit formulas for 
these will simply be long cumbersome analogs of those presented above. One question 
that can be asked is what are the generalized Lorentzian signature triangle inequalities 
for these Lorentzian simplices? One can show that the (4,1) simplex has its 4 time-like 
edge lengths completely unconstrained, provided its Euclidean tetrahedron is realizable, 
while the (3,2) simplex has its edge lengths subject to complicated constraints [34]. 
Thus, it is safe to say that the constrained Lorentzian configuration space is very 
different from that of the Euclidean configuration space in any dimension. 

8 Conclusions 

We have formulated a model of simplicial Quantum Gravity which is analgous to the 
Quantum Regge Calculus, but in the Lorentzian domain. This model is formulated 
by using an insight from the development of CDTs, which is to use the structure of 
Lorentzian manifolds to guide which geometries appear in the path integral. It is 
interesting to note that in CDTs this results in decreasing the set of allowed geometries 
to those which have a causal structure; however in our model we increase the number of 
geometries by modifying and relaxing the triangle inequalities. In 1+1 dimensions the 
removal of the triangle inequalities allows every dynamical edge length in the model to 
be integrated over in a completely unconstrained fashion. This allows more freedom 
for analytical calculation, which we have used to derive scaling relations for the pure 
gravity theory. Unfortunately obstacles remain for both full calculation of the path 
integral and also for calculation of loop amplitudes. In calculating loop amplitudes, 
the integrals that arise have no known closed form solution and block us from deriving 
the full result. Unless explicit results for these integrals are found, further analytical 
work will have to proceed in different directions. 
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Another aspect of this model is that, in contrast to CDT, the pure gravity path 
integral does not exhibit any critical behavior. This means that there is no way to 
take the usual sort of continuum limit of the model which could be compared to other 
approaches. This same aspect is found in the Euclidean version of the model, and this 
is resolved in that model by introducing higher derivative terms into the action [14]. 
The exact same procedure can be done in this model by using the angle formulas given 
in equations (2.9)-(2.11). This modified model will most likely need to be studied 
numerically through Monte Carlo simulations, as analytical work will run into the 
same obstacles found in the preceding sections. On the other hand, the lack of phase 
transition could be interpreted as a plus: A potential polymerized phase is excluded 
kinematically, while in 1+1 dimensions the occurrence of a rough/ fractal/ crumpled 
phase seems unlikely given the exact result (R) = and the well-behaved nature of the 
causal diamond truncation. This strongly suggests the theory is always in a smooth 
phase. Large universes (with spatial size controlled by the number N s of triangles in 
the spatial direction, and temporal size controlled by the number N t of triangles in the 
temporal direction) seem to be the norm rather than the exception. 

There are prospects for further analytical work using the causal diamond truncation 
introduced in section 4. Because in pure 1+1 gravity the only degree of freedom is the 
Liouville field, the causal diamond truncation could be justified since the varying space- 
like edge lengths act like a simplicial version of the Liouville field. Future work could 
focus on calculating quantities in this truncation such as the string susceptibility x<t>- 

The model is also safe from the occurrence of the spikes which trouble the Euclidean 
QRC In both the proper-time and the Causal Diamond truncation of the model the 
partition function with fixed space-time area is well behaved with the appearance of 
spikes completely suppressed. In the full model, there are regions of the configuration 
space in which the geometry (while not spiky) can become "pinched". These regions 
may or may not be important, however if they are it is not clear that they are patho- 
logical; instead they may correspond to an interesting phase of the model, perhaps only 
controllable by adding an R? term in the action. 

The model can be extended to higher dimensions, however generalized Lorentzian 
triangle inequalities are introduced, and difficulties analogous to the ones found in the 
Euclidean theory appear. The Lorentzian model in higher dimensions is not related 
to the Euclidean model by simple Wick rotation, as there are still sets of edge lengths 
which are unconstrained, and therefore the configuration space of the Lorentzian model 
is larger than that of the Euclidean one. The model in higher dimensions will need to 
investigated numerically, and it will be interesting to see how the results compare with 
those of CDT. 

In summary, in this paper we have demonstrated that Quantum Regge Calculus in 
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the Lorentzian domain differs from that in the Euclidean domain in ways considerably 
more complicated than by a simple Wick rotation. This is a specific case of the general 
fact that Lorentzian and Euclidean quantum gravity are not trivially connected. In 
order to study Lorentzian quantum gravity using QRC it is therefore necessary to take 
into consideration the differences we have outlined in this paper. 



References 

J. Ambjorn and R. Loll, "Nonperturbative Lorentzian quantum gravity, causality and 
topology change", Nucl. Phys. B 536, 407 (1998) [arXiv:hep-th/9805108]. 

R. Loll, "Discrete Lorentzian quantum gravity", Nucl. Phys. Proc. Suppl. 94 (2001) 96 
[arXiv:hep-th/0011194]. 

J. Ambjorn, A. Dasgupta, J. Jurkiewicz, R. Loll, "A Lorentzian cure for Euclidean 
troubles", Nucl. Phys. Proc. Suppl. 106 (2002) 977-979. [hep-th/0201104]. 

J. Ambjorn, J. Jurkiewicz, R. Loll, "Spectral dimension of the universe", Phys. Rev. 
Lett. 95 (2005) 171301. [hep-th/0505113]. 

J. Ambjorn, J. Jurkiewicz, R. Loll, "Reconstructing the universe", Phys. Rev. D72 
(2005) 064014. [hep-th/0505154]. 

J. Ambjorn, J. Jurkiewicz and R. Loll, "Quantum Gravity, or The Art of Building 
Spacetime" , arXiv:hep-th/0604212. 

F. David, "Simplicial quantum gravity and random lattices", arXiv:hep-th/9303127. 

D. Benedetti, J. Henson, "Imposing causality on a matrix model", Phys. Lett. B678, 
222-226 (2009). [arXiv:0812.4261 [hep-th]]. 

J. Ambjorn, J. Jurkiewicz and R. Loll, "The Universe from scratch", Contemp. Phys. 
47, 103 (2006) [arXiv:hep-th/0509010]. 

M. Visser, "Wormholes, Baby Universes, and Causality", Phys. Rev. D41 (1990) 1116. 

M. Rocek and R. M. Williams, "Quantum Regge Calculus", Phys. Lett. B 104, 31 
(1981). 

R. M. Williams, "Quantum Regge Calculus in the Lorentzian Domain and its 
Hamiltonian Formulation", Class. Quant. Grav. 3, 853 (1986). 

R. Sorkin, "Time Evolution Problem in Regge Calculus", Phys. Rev. D 12, 385 (1975) 
[Erratum-ibid. D 23, 565 (1981)]. 

H. W. Hamber and R. M. Williams, "Two-Dimensional Simplicial Quantum Gravity", 
Nucl. Phys. B 267, 482 (1986). 



- 34 - 



5] M. Gross and H. W. Hamber, "Critical properties of two-dimensional simplicial 
quantum gravity", Nucl. Phys. B 364, 703 (1991). 

3] H. W. Hamber, "Phases of simplicial quantum gravity in four-dimensions: Estimates 
for the critical exponents", Nucl. Phys. B400, 347-389 (1993). 

7} T. Regge, "General Relativity Without Coordinates", Nuovo Cim. 19, 558 (1961). 

3] R. M. Williams, "Recent progress in Regge calculus", Nucl. Phys. Proc. Suppl. 57, 
73-81 (1997). [gr-qc/9702006]. 

)] K. Tate and M. Visser, "Simplicial 1 + 1 gravity in terms of causal diamonds" , in 
preparation. 

]] H. W. Hamber and R. M. Williams, "On the measure in simplicial gravity", Phys. Rev. 
D 59, 064014 (1999) [arXiv:hep-th/9708019]. 

1] J. Ambjorn, J. L. Nielsen, J. Rolf and G. K. Savvidy, "Spikes in quantum Regge 
calculus," Class. Quant. Grav. 14, 3225 (1997) [arXiv:gr-qc/9704079]. 

2] E. Bittner, A. Hauke, H. Markum, J. Riedler, C. Holm and W. Janke, "Z(2) Regge 
versus standard Regge calculus in two-dimensions", Phys. Rev. D 59, 124018 (1999) 
[arXiv:hep-lat /9903028] . 

i] H. W. Hamber, R. M. Williams, "Higher Derivative Quantum Gravity On A Simplicial 
Lattice", Nucl. Phys. B248, 392 (1984). 

1] J. Ambjorn, J. Jurkiewicz, C. F. Kristjansen, "Quantum gravity, dynamical 
triangulations and higher derivative regularization" , Nucl. Phys. B393, 601-632 
(1993). [hep-th/9208032]. 

5] R. Nakayama, "2-D quantum gravity in the proper time gauge", Phys. Lett. B 325, 
347 (1994) [arXiv:hep-th/9312158]. 

3] D. S. Jones, "Incomplete Bessel Functions, I", Proc. Edinb. Math. Soc. 50, (2007). 

7} F. E. Harris, "Incomplete Bessel, generalized incomplete gamma, or leaky aquifer 
functions," J. Comput. Appl. Math 215, 1 (2008). 

3] X. Yang, "A matrix trace inequality", J. Math. Anal. Appl. 250, (2000). 

3] L. Liu, "A Trace Class Operator Inequality", J. Math. Anal. Appl. 328, (2007). 

]] J. Ambjorn, A. Gorlich, J. Jurkiewicz and R. Loll, "CDT - an Entropic Theory of 
Quantum Gravity," arXiv: 1007.2560 [hep-th]. 

[31] H. W. Hamber, "Simplicial Quantum Gravity", Les Houches Summer School 1984. 

[32] J. Ambjorn, J. Jurkiewicz and R. Loll, "Dynamically triangulating Lorentzian 
quantum gravity", Nucl. Phys. B 610, 347 (2001) [arXiv:hep-th/0105267]. 

[33] K. Wirth and A. S. Dreiding, "Edge lengths determining tetrahedrons", Elem. Math. 



- 35 - 



64 (2009) 160-170. 

[34] K. Tate, M. Visser, "Realizability of the Lorentzian (n,l)-Simplex," 
[arXiv:1110.5694 [gr-qc]]. 



- 36 - 



